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ABSTRACT 

The nucleus of the nearby galaxy, NGC 1097, is known to host a young, 
compact (r < 9 parsec) nuclear star cluster as well as a low-luminosity active 
galactic nucleus (AGN). It has been suggested both that the nuclear stellar cluster 
is associated with a dusty torus, and that low-luminosity AGN like NGC 1097 
do not have the torus predicted by the unified model of AGN. To investigate 
these contradictory possibilities we have acquired Gemini/T-ReCS 11.7 /im and 
18.3 fim images of the central few hundred parsecs of this galaxy at <45 parsec 
angular resolution, in which the nucleus and spectacular, kiloparsec-scale star- 
forming ring are detected in both bands. The small-scale mid-infrared (mid- 
IR) luminosity implies thermal emission from warm dust close to the central 
engine of this galaxy. Fitting of torus models shows that the observed mid-IR 
emission cannot be accounted for by dust heated by the central engine. Rather, 
the principal source heating the dust in this object is the nuclear star cluster 
itself, suggesting that the dust that we detect is not the torus of AGN unified 
schemes (although it is also possible that the dusty starburst itself could provide 
the obscuration invoked by the unified model). Comparison of Spitzer/IRS and 
Gemini/GNIRS spectra shows that, although polyaromatic hydrocarbon emission 
(PAH) bands are strong in the immediate circumnuclear region of the galaxy, 
PAH emission is weak or absent in the central 19 parsecs. The lack of PAH 
emission can probably be explained largely by destruction/ionization of PAH 
molecules by hard photons from the nuclear star cluster. If NGC 1097 is typical, 
PAH emission bands may not be a useful tool with which to find very compact 
nuclear starbursts even in low-luminosity AGN. 

Subject headings: dust, extinction - - galaxies: active - - galaxies: nuclei - 
infrared: galaxies — galaxies: individual (NGC 1097) 

1. Introduction 

The dusty molecular torus, cornerstone of the unified model of active galactic nuclei 
(AGN), has been the subject of much observational and theoretical scrutiny in recent years. 
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The detection of broad emission lines in polarised light in a number of type 2 objects 


(Antonucci & Miller 


1985: 


Miller & Goodrich 


1990: 


Young et al. 


1996: 


Moran et al. 


2000; 


Iran 


2001 


) demonstrates the overall validity of the unified model, and structures thought to 



correspond to t he torus itself hav e now been resolved through mid-IR interferometry of one 



Seyfert galaxy (IJaffe et al.l 12004 ). However, questions remain as to the origin and precise 
nature of the obscuration and the extent to which it can account for the differences between 
the many classes of AGN. 

In particular, the possibility that the torus and star formation are intimately connected 
has been raised by several authors on theoretical grounds. For instance, the material in the 
nucleus needs to have a certain scale height to obscure the nucleus and account for the relative 
numbers of types 1 and 2 AGN, and it has been suggested that the energy necessary to main- 
tain this thickness might be supplied by supernovae in nuclear starbursts (IWada fc Norman 
20021 ). Indeed, some or all of the nuc lear obscuration has been proposed to arise in dusty 



clouds ejected from such starbursts (jCid Fernandes fc Terlevic 



1995 



Fabian et al. 1998 



Watabe fc Umemural 120051 ) or in optically thick stellar winds (IZier fc Biermannl |2002| ) . Ob- 
servationally, young stars are found quite commonly in the inner few hundred pc of AGN (e.g. 



Gonzalez Delgado et alll998l : IStorchi-Bergmann et all2000l : Ilmanishill2002l . 120031 ; iRodrfguez-Ardila fc Viee 
2OO3I). 



However, while few-hundred-parsec-scale star formation is often observed in AGN, few 
young stellar clusters have to date been detected at smaller scales. The outer diameter of 



resolved in single-dish 10-20 /zm observations of even the nearest objects (e.g. 


Bock et al. 


2000; 


Radomski et al. 


2003; 


Soifer et al. 


2003; 


Packham et al. 


2005; 


Mason et al. 


2006), and 


interferometric measurements of NGC 1068 


suggest that its mid-IR-emitting regions are 



only -3 pc (~0.04' 
the Circinus galaxy ( 



in diameter (Jaffe et al.l I2QQ4T ) . Near-IR adaptive optics imaging of 



Prieto et al.ll2004h resolves a source of FWHM~2 pc, consistent with a 



parsec-scale torus at those wavelengths. 



The discovery by IStorchi-Bergmann et al.l (120051 . hereafter SB05) of UV absorption sig- 
natures of a 1O 6 M , 10 6 year-old stellar cluster within only 9 pc of the nucleus of the nearby 
LINER/Seyfert 1 galaxy, NGC 109^, is therefore of obvious interest. The star cluster — but 
not the broad emission lines from the nucleus itself — suffers about 3 magnitudes of visual 
extinction, consistent with the possibility of starburst activity physically associated with a 



1 Although NGC 1097 hosts a classic LINER nucleus ([Phillips et al.lll984l ). broad , double-peaked emission 
lines in its spectrum are characteristic of rotating gas in an AGN accretion disk (jStorchi-Bergmann et al. 
Il993h . 
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dusty AGN torus. Young superstellar clusters or recent st ar formation are also known to 



exist in the central few parsecs of a handful of other AGN (IColina et al.l 120021 ; iDavies et al. 
2006bl : iMueller Sanchez et al.ll2006t ). 



While the results of SB05 are suggestive of a starburst genuinely associated with an 
obscuring torus, the observation that many LINERs host nuclear UV p oint sources raise s 
some doubt as to whether LINER- type nuclei in fact host obscuring tori (IMaoz et al.ll2005l ). 
This is expected on theoretical grounds in the context of the disk-wind scenario for the torus: 
below a certain luminosity the accretion on to the black hole can no lo nger sustain the mass 
outflow necessary to obscure the nucleus (jElitzur fc Shlosmanl 120061 ). Moreover, how the 
properties of LINER tori, if they exist, might relate to those observed in Seyfert galaxies is 
completely open to question. Further observations of LINER nuclei to probe their detailed 
nuclear structure are therefore desirable. 

In addition, the UV evidence of star formation detected by SB05 is very susceptible to 
extinction, so UV detections of starbursts in AGN tori are likely to be limited to cases of 
nuclei in which the stellar clusters are not too deeply buried within the nuclear dust. Inter- 
preting the energetics of UV-detected starbursts will be complicated by uncertain extinction 
corrections. A more useful tracer of the nuclear starburst phenomenon would be an infrared 
star-formation indicator that is both unaffected by dust extinction and accessible from the 
ground, to take advantage of the spatial resolution achievable from ground-based telescopes. 
The 3.3 /im polycyclic aromatic hydrocarbon (PAH) emission band fits both of these crite- 
ria. Usually attributed to vibrational relaxation of aromatic hydrocarbon molecules excited 



in extraga 


— „„, ~ 
actic environments (e.g. 


Moorwood 1986: Roche c 


t al. 


1991: 


Genzel et al. 1998: 


Lutz et al. 


1998; Rieopoulou et al. 


1999; 


Laurent et al. 


2000; 


Trar 


i 2001; 


Imanishi & Wada 


2004: Peeters et al. 2004 : iRisaliti et al. 


2006). 



The suite of PAH emission features is commonly detected in large-a perture ISO spec- 



tra o f nearby active galaxies with circumnuclear star-forming regions (IRigopoulou et al. 
19991 ). but on moving to apertures of 3-5" (~ few hundred p c) or less, the bands are 



often weak or abse nt (jRoche et al. 1991 ; Le Floc'h et al. 2001 ; Siebenmorgen et al. 2004 
Roche et all l2006al fbh . Modelling suggests that the carriers of th e bands are easily de- 



stroyed by the harsh radiation field c lose to the active nucleus itself (jRoche fc Aitken 



Siebenmorgen et al 



For instance, 



12004 



1985 



Voit 



992). However, this appears not always to be the case. 
Imanishi fc Wadal (120041 ) find that the surface brightness of the nuclear PAH 
emission in Seyfert 1 galaxies is comparable to that in starburst galaxies, and the 3.3 /im 
PAH band has been shown to be present in several Seyfert galaxies of both types 1 and 
2 on fairly small spatial scales, including, in the case of NGC 3227, within only 60 pc of 
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the central engine (IRodriguez-Ardila fc Viegaal2003l ). To exist at such small distances from 
these AGN, it seems likely that the fragile PAH band carriers are shielded from the nuclear 
X-rays. This intriguing result could perhaps be explained if the carriers exist within the 
dense, dusty clouds of the torus itself. 

To further investigate this possibility, we have performed L band (~3-4 jttm) spec- 
troscopy of the nucleus of NGC 1097 (D=17 Mpc for H = 75 km sec" 1 Mpc" 1 ; 1"= 85 pc) 
to search for the 3.3 /im PAH emission band. We have also obtained images at 11.7 and 
18.5 /im, to search for the compact mid-IR source expected to be a signature of the torus 
and to constrain its properties. The data were taken using the Thermal Region Camera 
and Spectrograph (T-ReCS) and the Gemini Near-Infrared Spectrograph (GNIRS) on the 
Gemini South 8.1 m telescope, yielding <0.52" spatial resolution and therefore probing dis- 
tances < 22 pc from the nucleus. We also discuss complementary Spitzer Space Telescope 
Infrared Spectrograph (IRS) observations of the central 3.6 x 3.6" (approx. 300 x 300 pc) 
of this galaxy. $2] describes all of these observations and the reduction of the data. In $3] we 
present the nuclear data and describe the characteristics of the compac t source dete c ted in 
the mid-IR, and compare our su barcsecond-reso l ution data and that of IPrieto et al.l (120051 ) 
with the clumpy torus models of iNenkova et al.l (120021 ). The conclusions that can be drawn 
from this work are discussed in $H 



2. Observations and data reduction 
2.1. Mid-IR imaging 

Narrow-band N (11.7 /im) and Q (18.3 /im) band images of NGC 1097 were taken on 



2005 September 16 and 17 using the Gemini South mid-IR instrument T-ReCS ( jTelesco et al. 



19981 ). T-ReCS uses a Raytheon 320 x 240 pixel Si: As IBC array, providing a pl ate scale of 



0.089 " per pixel. The detector was used in correlated quadruple sampling mode (jSako et al. 



20031 ). Images were obtained in the 11.7/im (5X = 1.13/im, 50% cut-on/off) and 18.3/im (5\ 
= 1.51/im, 50% cut-on/off) filters using the standard chop-nod technique to remove time- 
variable sky background, telescope thermal emission and "1/F noise. The chop throw was 
15" and the telescope was nodded 15" in the direction of the chop every 30 s. T he chop angle 
was chosen based on the 5"-aperture N band data of iTelesco fc Gatleyl (Il98ll ) to minimise 
the extended emission in the off-source beam and was fixed at 90 ° east of north. The total 
on-source time for the 11.7 /im and 18.3 /tm observations was 456 s and 912 s respectively 

The data were reduced using the Gemini IRAF package. The difference for each chopped 
pair for each nod-set was calculated, and the nod-sets then differenced and combined until a 



- 6- 



single image was created. Chopped pairs obviously compromised by variable sky background, 
high electronic noise, or other problems were searched for but no data needed to be discarded. 

The Cohen standards SAO 216405 and SAO 193679, and PPM 278407 were observed 
for flux and PSF calibration. Observations of the standards taken approximately 1 hr apart 
showed a variation in counts of around 10%, typical of mid-IR photometry. The PSF ob- 
servations were made immediately before or after the NGC 1097 observations and with the 
same instrumental configuration. The FWHM of the PSF star was 0.41" (standard de- 
viation ~0.03") at 11.7 /im and 0.52" in the 18.3 fim filter. Short PSF or flux standard 
observations are comparable to longer source observations as the closed-loop active optics of 
Gemini provides a similar PSF when taken at a similar telescope pointing and time (see e.g. 
Radomski et al.ll2003l ; iPackham et al.ll2005l although seeing variations may of course affect 
the final observed PSF). As judged by the FWHM of the standard stars, the PSF was stable 
over the course of these observations. 



2.2. L-band spectroscopy 

L band spectroscopy of the nucleus of NGC 1097 was obtained using GNIRS, the near-IR 
spectrograph on the Gemini South telescope, on 2005 November 13. The short red camera 
and 3-pixel (0.45") slit were used together with the 32 1/mm grating, resulting in R~1100. 
The total integration time was 30 minutes, and the telescope was nodded by 20" every 30 
seconds in an ABBA pattern so that the sky background emission could be subtracted out. 
Sky conditions were photometric during these observations. 

The nucleus was acquired in the H band to take advantage of the contrast between 
object and sky at those wavelengths. A point source contributes 52% of the H band flux in a 



0.4" diameter aperture, that fraction increases with wavelength (jPrieto et al.ll2005l ). and we 
observe an unresolved source at 11.7 fim (§ I3.ip . suggesting that the location at which the 
slit was centered at 1.6 /im also corresponds to the flux peaks at 3.3 and 11.7 /im. The slit 
was orientated at 9 ° east of north in order to encompass one of the circumnuclear emission 
knots visible in the T-ReCS image, but that knot proved too faint to yield a useful spectrum 
and will not be considered further. 

The data were reduced using the Gemini IRAF package in combination with standard 
IRAF and Starlink FIGARO routines. After bad pixel masking, flatfielding and subtracting 
in pairs, the individual sky-subtracted spectra were averaged together, with the positive and 
negative beams cross-correlated and shifted in both the spatial and dispersion directions 
before combining. The central three pixels (0.45") of the galaxy spectrum were extracted 
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using the standard star spectrum to define the curvature of the extraction aperture across 
the array. 

Wavelength calibration was accomplished using the numerous telluric absorption lines 
in the spectrum. The galaxy spectrum was then cross-correlated with the spectrum of 
the standard star, BS 466 (F7V), and divided by that spectrum. In order to improve the 
cancellation of telluric lines, the resulting spectrum was multiplied by a logarithmically- 
scaled model atmospheric spectrum^], binned to the appropriate spectral resolution. This 
had the effect of slightly reducing the amplitude of some of the residuals in the 3.3 /im 
region. The spectrum was then multiplied by a 6240 K blackbody curve. Flux calibration 
was with reference to BS 466, using a V band magnitude of 6.24 and V-L = 1.36. No attempt 
was made to correct for uncertain differential slit losses between the galaxy and standard 
star (due to e.g. seeing variations and different intrinsic spatial flux distributions) and, as 
is usual for flux calibration of narrow-slit spectral data, the uncertainty of the derived flux 
density is considerable. 



2.3. IRS spectroscopy 

Publically-available spectra of the nucleus of NGC 1097 taken with the short-low (first 
order) module of the Spitzer Space Telescope's Infrared Spectrograph (IHouck et al.l 12004 
IRS) were retrieved from the Spitzer data archive (PID 14; PI Houck). These R~100 ob- 
servations represent a total integration time of 56 s. Files output from the IRS pipeline at 
the Spitzer Science Center (which includes ramp fitting, dark-sky subtraction, droop cor- 
rection, linearity correction, flatfielding and wavelength calibration) had background light 
subtracted by subtracting the first-order slit exposed during second-order observations of 
NGC 1097 taken a few minutes previously. Spectra were extracted from the background- 
subtracted, coadded frames using Spice, the SPItzer Custom Extractor. The extraction 
window scales with wavelength but was defined to be 2 pixels (3.6") at 11.2 /im, equal to 
the 3.6" slit width and comparable to Spitzer's diffraction limit at 11.2 /im (3.3"). 

The spatial profile of the IRS spectrum shows flux between the nucleus and the circum- 
nuclear ring ~10" distant, presumably from the wings of the PSFs of these bright sources 
combined with real diffuse emission from that region. Neither the point source nor the ex- 
tended source flux calibration provided by Spice is optimised to deal with this situation. 
However, the conclusions drawn from these observations do not depend critically on a highly 
accurate calibration, so we simply normalise the mean flux in the IRS spectrum to that given 



Lordl (|l992h . obtained from www.gemini.edu/sciops/ObsProcess/obsConstraints/ocTransSpectra.html 
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by lTelesco fc Gatleyl (ll98ll ) over a similar bandpass in a slightly larger aperture (5" diameter 
vs. 3.6 x 3.6"). 



Results 



3.1. The mid-IR images and model fits 

The T-ReCS mid-IR images of NGC 1097 are shown in Figures [T] and HJ A compact 
nuclear source and the well-known circumnuclear star-forming ring are clearly detected at 
both 11.7 and 18.3 /mi. The starburst ring emits strongly in the mid-IR, with each of the 
three knots in the northern part of the ring contributing a similar level of flux to that of 
the nucleus itself. The 15" maximum chop throw of the Gemini telescopes means that the 
positive and negative images partially overlap in the region of the circumnuclear ring. Inside 
the ring, faint extended emission is very probably present but is difficult to detect without a 
clear sample of the sky on which to measure the zero level of the background. The flux of the 
nuclear point source is approximately 48 mJy at 11.7 /im and 64 mJy at 1 8.3 /im (Tabl et!]), 
rough ly 75% and 27% of the broad-band, 5"-aperture measurements of iTelesco fc Gatley 
( 1l98ll ) at N and Q respectiv ely, and consistent with the 60 mJy, 1.5"-aperture N-band upper 
limit of lGorjian et al.l (120041 ) . Faint emission between the nucleus and ring likely accounts for 
the difference between our data and that of ITelesco fc Gatleyl (Il98ll ). most clearly affecting 
the Q band photometry. 

Table [1] gives the flux densities of the three regions of emission in the circumnuclear 
ring that are clearly detected at 18.3 /im and least likely to be affected by the negative chop 
beams. The aperture sizes used were determined partly by the proximity of neighboring 
emission regions and/or the edge of the detector. The N/Q flux density ratios are in the 
region of 0.4 for the ring, but rather higher, about 0.75, for the nucleus. 

The FWHM of the central compact source was measured to be 0.43" at 11.7 /im, com- 
parable to the mean FWHM of the comparison star observations (0.41±0.03", § [2]). This 
translates to a dia meter of < 37 pc, consistent with the NIR limit of 10 pc diameter of 
Prieto et al.l (120051 ) . Thus, like all AGN observed non-interferometrically to date, the nu- 
clear source of NGC 1097 is not resolved at its FWHM in the N band. Limited S/N prevented 
a meaningful measurement of the source size at 18.3 /im. 
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Table 1. Photometry of the nucleus and circumnuclear ring of NGC 1097 



Location 1 


Aperture diameter 
(arcsec) 


11.7 (im flux density 
(mjy) 


18.3 fim flux density 
(mjy) 


Nucleus 


3.0 2 


48 


64 


A 


1.6 


34 


83 


B 


2.0 


75 


169 


C 


1.6 


25 


63 



!See Fig. [2] 

2 An aperture of 3" diameter was used for the photometry of the nucleus, chosen 
based on the profiles of the standard stars to ensure that all of the flux from the 
unresolved source was collected. A surrounding "sky" annulus was used to subtract 
off any residual sky background and diffuse emission surrounding the nucleus. The 
(unknown) spatial variation of surrounding diffuse emission will therefore determine 
the degree to which the photometry represents the unresolved source alone, but as 
any such emission must be close to the sensitivity limit of the instrument, we expect 
it to make a negligible contribution to the measured nuclear flux. 
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Fig. 1. — Central portion of the T-ReCS 11.7 image of NGC1097, smoothed with a 2-pixel 
Gaussian. The (unguided) negative off-source images are evident at either side of the positive 
galaxy image. 
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Fig. 2. — As for Figure 1, but at 18.3 /zm. Photometry of clumps A, B and C is given in 
Table EJ 
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While AGN- or star-formation-heated dust is expected to emit in the mid-IR, other 
processes can a lso c ontrib ute to the nuclear mid - IR em ission. For instance, as shown by 
Perlman et al.l (120011 ) and IWhysong fc Antonuccil (12004 ). the 10 /im emission from M87 on 
subarcsecond scales can be fully accounted for by synchrotron radiation from the inner regions 
of the ra dio jet. However , the t hin disk + jet + radiatively inefficient accretion flow (RIAF) 
model of iNemmen et al.l (120061 ). which does not include any thermal dust emission, predicts 
mid-IR luminosity around an order of magnitude lower in both the N and Q bands than that 
observed from the unresolved nuclear source in NGC 1097 (log[z/L^(erg s -1 )] ~ 41.6 at both 
N and Q for D=17 Mpc. Altering the mass accretion rate or other parameters of the thin 
disk + RIAF model cannot account for the observed mid-IR flux, given the observational 
constraints at shorter wavelengths. Even allowing for the fact that the jet parameters in the 
Nemmen et al. model are not well constrained, we conclude that the bulk of the mid-IR 
emission that we detect comes from warm dust in the region of the central engine. 

To investigate the contribution of the AGN to the heating of the nuclear dust, we fit 
the observed fluxes with the inhomogeneous torus models of Nenkova et al. (2002). In these 
computations, individual clouds are optically thick, and the inner radius of their toroidal 
distribution is located at the dust sublimation radius. While the total extent of the torus 
may reach 30 to 100 times the inner radius, the smallest physical scales, where the clouds are 
hottest, dominate the emission at near/mid-IR wavelengths. The physical parameters that 
describe the distribution of clumps include the average number of clouds along radial rays 
in the equatorial plane, N , and the optical depth per cloud in the V band, Ty. The clouds 
follow a power law distribution in radius, r, with number declining as r~ q . The mid-IR 
emission is not sensitive to the outer radius of the cloud distribution, so we fix it at 30 times 
the inner radius. 

We require that the clump distribution give a high probability (> 0.5) of direct views of 
the AGN, given the l ack of reddening of the broad Ha line and the l ow X-ray column density 



towar d the nucleus (ITerashima et al.l |2002| ; INemmen et al.l 120061 ) . IStorchi-Bergmann et al. 
( 2003 ) find the inclination angle of the accretion disk i = 34^2 . We assume that the torus is 
aligned with the accretion disk and adopt i = 30° in the coarse model grid. The torus does 
not have a sharp edge, but rather a Gaussian angular distribution away from the equatorial 
plane of width a. In terms of the angle from the equatorial plane, (3, the number of clouds 
along a ray N(j3) = iVoexp(— /3 2 /cr 2 ). To allow direct views to the nucleus with the viewing 
angle constrained, we consider a = 30°, which roughly corresponds to a torus half-opening 



3 Note that contrary to the description in their paper, Storchi-Bergmann et al. ( 2003h intend i to have its 
usual meaning, measured with respect to the rotation axis (not the plane) of the disk, so edge-on views have 
i — 90° (M. Eraclcous, private communication). 
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angle of 60°, and a = 45° (for half-opening angle ~ 45°). With the broader distribution, the 
probability of viewing the unobscured nucleus is generally low, and only N < 2 is acceptable 
at i = 30°. 

We consider three different heating spectra and fit the sum of the intrinsic spectrum and 
the resulting dusty emission to the ratio of our N- and Q-band fluxe s. We scale the results 
to the flux at Q (Figs. [3] and H}, except when using the result of iNemmen et al.l (120061 ) 
for the input spectrum, in which case multi-wavelength observations already constrain the 
luminosity. Also shown in Figures and H] are "AGN magnitudes," the J, H and K band 
flux densities of the nuclear point source derived by fitting unresolved emission and an 
underlying galaxy component to near-IR imaging data ( Prieto et al. 2005 1 4 ! We speculate 
that NGC 1097 will remain unresolved in the mid-IR at resolutions similar to the PSF 



used by iPrieto et al.l (120051 ) to fit the near-IR data. Such a result would imply that the 
mid-IR emitting r egion is 5-10 pc i n radius, as has been found for other nea rby AGN (e.g. 
Soifer et al.l 120031 : Ijaffe et all 12004! : IPackham et~all I2OO5I : iRoche et al.l l2006bh . We therefore 
believe it likely that these near-IR data sample the same spatial region as our mid-IR points. 

The first dust-heating spectrum we consider is that of a standard AGN (in which \F\ oc 
A~ ' 5 from 0.1 to 1/rni). Although the low observed optical and UV fluxes together with 
the minimal obscuration along the line of sight rule out this model for NGC 1097, it serves 
as a useful point of comparison. In this case, the dusty material dominates the mid-IR 
emission, even when the AGN is viewed directly, and the AGN itself dominates at optical 
and near-IR wavelengths. The dust sublimation radius scales with the square root of the 
bolometric luminosity of the heating source and is a function of the heating spectrum. This 
AGN heating spectrum and the observed = 10 42 erg s _1 set the inner radius at 0.01 
pc. The observed N/Q ratio and a likely unimpeded view of the central engine exclude 
q = 1, 2, and 3 radial distributions. Two reasonable models both have q = 0, o = 30°, and 
iVo = 2, with t v = 40 and 100 in models "AGN 1" and "AGN 2," respectively. Scaling 
to the observed flux at Q yields Lboi = 4 x 10 42 erg s _1 , which is a few times higher than 
the observed AGN luminosity. More severely, they significantly exceed the observed optical 
emission, even after correcting for the measured extinction. This implies that although dust 
heating by a standard AGN radiation field can reproduce the mid-IR SED of NGC 1097 
well, this cannot in fact be the dominant dust-heating mechanism in this galaxy. 



4 The values given in IPrieto et al 
(private communication) . 



(|2005l ) are erroneous; corrected values were supplied by M. A. Prieto 
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Fig. 3. — Inhomogeneous dusty torus models, fit only to the 11.7 and 18.3 /im fluxes (squares) 
and constrained to allow a direct view of the central engine. The sum of the intrinsic 
heating spectrum and the resulting dust emission is plotted as a solid line in each case. 
Several different obscuring geometries could produce the mid-IR emission if a standard AGN 
spectrum (black and dark gray lines) heats the dusty clouds. These AGN models are scaled 
to the 18.3 /im flux, and applying foreground extinction of the intrinsic sour ce of Ay = 0.12 



mag r esults in the corresponding dotted lines. Alternatively, using the fit of iNemmen et al. 



(120061 ) as the heating spectrum, none of the allowed torus geometries gives a good fit to the 
mid-IR flux ratio or the observed flux in either band. We plot the b e st of these (light gray 
line), fixed to the intrinsic luminosity determined by INemmen et al.l (120061 ). None of these 
AGN or RIAF solutions remains consistent with the observations at all wavelengths. We 
plot a smoothed version of the L-band spectrum (thick line) and small-scale near-IR 'AGN 
magnitudes" (inverted triangles; see text). We also show the total observed emission (filled 
triangle) and reddened starburst contribution (open triangle) that SB05 measure at 9000A. 
ISO and IRAS data, which could in principle constrain the models in the far-IR, are not 
shown as the large apertures involved contain a very large contribution from emission from 
the bright, kiloparsec-scale starforming ring. See the text for a complete description of the 
different models. 
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Nemmen et al.l (120061 ) present the broad-band SED of the nuclear region of NGC 1097. 
The data reveal the intrinsic spectrum at optical and UV wavelengths, which is responsible 
for heating the surrounding dust, although they do not constrain it well in the IR. Motivated 
by the empirical result at short wavelengths, we consider Nemmen et al.'s thin disk + jet + 
RIAF model fit to the data as an alternate heating source of the nuclear dust. We use this 
spectral model (which we refer to as the "Nemmen06" model for simplicity) as the input to 
the clumpy torus distribution of clouds. 

The Nemmen06 spectrum is weaker than the standard AGN at wavelengths below lfim 
and relatively stronger than the standard AGN at longer wavelengths. The UV-weak con- 
tinuum heats the clouds inefficiently, so the dust sublimation temperature is reached closer 
to the central source than in the standard AGN model of the same luminosity. With the 
Nemmen06 heating and = 10 42 erg s -1 , the inner radius is located at 0.003 pc, while 
still scaling with the square root of luminosity. Fixing the viewing angle (i = 30°), none of 
the a = 30° models fits the data. Several of the a = 45° models with q = or 1 agree with 
the mid-IR flux ratio. However, the required direct view to the AGN rules out these models, 
which have Nq = 20 so the probability of an unobscured view of the nucleus is only 0.03. For 
comparison with the standard AGN results, we show the best of the allowed models based 
on the Nemmen06 spectrum, in which q = 0, a = 45°, iVo = 2, and ry = 60. We adopt 
the observed luminosity of the input Nemmen06 spectrum and therefore do not scale the 
result to the mid-IR data as we did with the standard AGN model. The resulting model, 
including the dust emission, underestimates the flux density at Q by a factor of two, while 
overestimating the flux density at L by a factor of two. 

Both the standard AGN and Nemmen06 spectra produce similar SEDs of the torus 
alone, with a peak in F\ around 20-40/iin. However, because the total model is the sum of 
the input and torus spectra, the significant differences between the input spectra strongly 
affect the relative mid-IR contribution of the torus. The intrinsic AGN spectrum is very 
weak in the mid-IR, so the total mid-IR emission of this model is essentially that of the 
torus alone. Here the torus is thousands of times brighter than the input spectrum in the 
mid-IR, even when the torus scale height is small. In contrast, the Nemmen06 continuum is 
strong at near- and mid-IR wavelengths, and its torus does not provide more than 3 times the 
flux input at 11.7 or 18.3/im. Consequently, these models generally overpredict the observed 
N/Q flux ratio. In the Nemmen06 case, the entire resulting optical to mid-IR spectrum is 
very sensitive to the input spectrum. The original model was not constrained well in the IR, 
and even without dust reprocessing it overpredicts the observed near-IR nuclear emission. 
A more severe problem of the Nemmen06-heated models is that the optical/mid-IR flux 
ratio is always too high. The optical/mid-IR signature of the intrinsic spectrum persists 
here, even when the dust-reprocessed contribution is included. Fundamentally, because the 
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short-wavelength flux is relatively weak, it cannot produce the strong mid-IR emission we 
measure. Thus, we conclude that the complete Nemmen06 model does not describe the IR 
data well. Changing the input contribution of the RIAF spectrum alone and neglecting dust- 
reprocessed emission is not a viable solution. The family of RIAF models cannot produce the 
observed mid-IR flux density, even when the underlying physical parameters of the models 
are varied. 

That the clumpy torus models with neither a standard AGN spectrum nor the Nem- 
men06 spectrum are able to simultaneously reproduce the mid-IR and the near-IR/optical 
emission of NGC 1097 suggests that the nuclear starburst contributes to the intrinsic lumi- 
nosity that the dusty clouds reprocess. In fact, a 10 6 MfD, 1 6 year-old starburst scaled by 
0.85 (SB05) has Lj, i = 5 x 10 42 erg s" - 1 jLeitherer et al.lll999h . several times higher than the 



AGN. To test this possibility we have adapted the clumpy torus models to use the young 
starburst spectrum for their input source, as an approximation to the mixed regions of dust 
and young stars presumably present in the nucleus. 

We use a spherically-symmetric distribution of clouds with N fixed at 0.7, for P> 0.5 of 
a direct view of the central engine. The observed N/Q ratio rules out q = 2 and 3 solutions, 
but we find good solutions with q — - 1 and ry = 20 - 60: q = and Ty = 20 ("SB 1" 
model), q = and t v = 40 ("SB 2" model), and q = 1 and t v = 60 ("SB 3" model). These 
three models have L = 1.2, 1.2, and 1.6 x 10 42 erg s _1 , respectively, when the model output 
is scaled to the observed mid-IR flux. This is somewhat lower than the starburst luminosity 
of SB05 but we note that that value relies heavily on the large extinction correction of the 
UV spectrum. Because we fit only the mid-IR data, changing the extinction would not alter 
our preferred model fits, although it would change the modeled observable optical emission 
(dotted lines in Figure HJ). The starburst heats the dust effectively because of its high 
UV luminosity, but its weakness in the infrared relative to both the standard AGN power- 
law and RIAF spectra means that the near-IR emission is predominantly from the dusty 
clouds, rather than from a bright AGN component that overpredicts the emission at optical 
wavelengths. The geometry of the starburst models is somewhat artificial, and we have no 
direct measurement of the spatial relationship between the starburst-heated clouds and the 
accretion disk. All the flux measurements include both the starburst and the central engine 
along with the emission they reprocess on scales of 10 pc. Nevertheless, basic considerations 
of energetics and the broad-band spectra of heating sources, as well as the specific model fits 
suggest that the nuclear starburst is the dominant contributor to the heating of the dust in 
the central few tens of parsecs of the galaxy. 

More detailed modelling would show more precisely the roles of the starburst and AGN 
in the heating of the dust. For example, we have opted against including any contribution 
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heated by the central engine while fitting the starburst to avoid the near-IR problems inherent 
in the underlying RIAF spectrum. In this simplified model fitting, we have used only the N/Q 
ratio to determine the model parameters. A further complication that is beyond the scope 
of this work would be to fit a broader spectral range and to include luminosity contraints 
simultaneously. Nonetheless, even the N/Q ratio by itself robustly excludes large portions 
of parameter space. We emphasize that whatever its ultimate power source, the mid-IR 
emission alone demonstrates that warm dust is present in the nucleus of NGC 1097. The 
dust is local to the nucleus and must obscure it from certain lines of sight whilst leaving the 
broad lines unobscured from our viewing angle, fulfilling at least two of the requirements of 
the AGN unified scheme. 



-18- 




SB 1 

SB 2 
SB 3 



Wavelength (|J.m) 



Fig. 4. — Fits of the starburst-heated dust model to the SED of NGC 1097. Solid lines show 
the sum of the intrinsic spectrum and the dust-reprocessed emission. Dotted lines show 
extinction of the starburst spectrum by Ay = 3. The starburst emission peaks in the UV, 
so it heats the dust effectively while not contributing directly to the net IR flux. Data as for 
Fig. El and model details are listed in the text. 
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3.2. The spectra 

The IRS short-low spectrum of NGC 1097 (Fig. [5]) shows a number of prominent emission 
features, including the 7.7, 8.6, 11.2 and 12.7 /zm PAH bands (this last line blended with 
the 12.82 fim [Ne II] fine structure line). NGC 1097 is therefore one of the growing number 
of AGN in which strong PAH emission is detected within as little as 150 pc of the active 
nucleus. 
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Fig. 5. — Spitzer IRS spectrum of the nucleus NGC 1097. Lines are drawn at 7.7 jum (PAH), 
8.03 /im (H 2 0-0(S4)), 8.6 /im (PAH), 9.66 /mi (H 2 0-0(S3)), 11.2 /mi (PAH), 12.28 /iin (H 2 
0-0(S2)),12.7 /im (PAH) and 12.82 /i m f[Ne III). The peak wavelen gths of the PAH bands 
vary somewhat from source to source (Ivan Diedenhoven et ajj 12004 ). 
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To estimate the 3.3 /zm band flux that would be expected from the 3.6 x 3.6" (approx . 
300 x 300 pc) region covered by the IRS spectrum, we use the finding of iHony et al.l (120011 ) 
that the 11.2 /im band is typically a factor of 3-4 stronger than the 3.3 /iin feature in a 
wide range of environments. The flux in the 11.2 /im band, 7.1 x 10 _13 erg sec _1 cm -2 (as 
obtained by a spline fit to the continuum and direct integration under the band), then 
implies ~ 2 x 10~ 13 erg sec _1 cm -2 in the 3.3 jum feature. The corresponding luminosities 
are 2.5 x 10 40 erg sec -1 and ~ 7 x 10 39 erg sec -1 for the 11.2 and 3.3 //m bands, respectively, 
within about a factor of two of the 3.3 /im PAH lumin osity of Seyfert 1 galax ies such as 
NGC 4235 and MCG -2-33-34 on similar physical scales Jlmanishi fc Wadall2004h . 

In the GNIRS spectrum of the central 0.45 x 0.45", however, no 3.3 jum feature is 
evident (Fig. Using the 95% confidence interval on the difference between the means 
of points near the peak of and just outside the featu re (and assuming a Gaussian profile 
with FWHM=0.04 /im; Ivan Diedenhoven et al.l 120041 ). we estimate the 3.3 /im PAH flux 
to be < 9 x 10 _15 erg sec -1 cm -2 . This is well below the flux that would be expected if 
the carriers of the 11.2 /im band were concentrated in the unresolved mid-IR source, as 
might be expected if they existed predominantly in warm, dense, dusty nuclear clouds, even 
allowing for uncertainties in the flux calibration of the GNIRS and IRS spectra. The spatial 
distribution of the mid-IR and PAH emission is very different. 
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Fig. 6. — L-band spectrum of NGC 1097. The black line represents the strongest 3.3 /im 
PAH band likely to be present in the spectrum (see text). 
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Comparison of the N band an d L' band continuum fluxes in similar apertures (ITelesco fc Gatley 
198ll : lAlonso-Herrero et alJlioOll l , combined with the above 11.2/3.3 /im band ratio, suggests 



that 3.3 /im band emission would peak at nearly twice the continuum level in our spectrum 
if the PAH band- and continuum-emitting material were similarly distributed over the cen- 
tral region of NGC 1097. That would imply an equivalent width (EW3.3) ~40 nm, whereas 
the maximum likely 3.3 /im feature flux estimated above implies EW3.3 <7 nm. This puts 
the EW3.3 of the nuclear starburst in NGC 1097 at th e low end of the values reported for 
few-hun dred-parsec-scale s t arbur sts by llmanishil ( 120031 ). although comparable to several of 
those in llmanishi fc Wadal (120041 ) . The lack of PAH emission in the nucleus implies an in- 
crease in the relative level of featureless L-band continuum in this aperture, possibly due to 
destruction of the PAH molecules (and/or dehydrogenation and ionization, whi ch affect the 



mole c ules' ability to emit at 3 .3 /xm and their likelihood of photodissociation; lAllain et al. 



19961 ; lAUamandola et al.lll999l ). This is examined in more detail in 



Discussion 



4.1. The nuclear mid-IR emission 



The detection of UV va riability, hard X-ray spectra and polarized broad emission lines 



in a number of LINERs (e.g. Barth et al.lll999l : iTerashima fc Wilson! 120031 : Maoz et al.ll2005f ) 



has demonstrated that at least some of these objects harbor accreting supermassive black 
holes at their centers. An obvious question, then, is whether these LINERs possess the dusty 
tori postulated by the AGN unified scheme. 

In the case of NGC 1097, the detection of unresolved thermal emission implies that 
there is warm dust in the nucleus of this galaxy that would obscure the nuclear broad lines 
if viewed from certain angles. As with Seyfert nuclei, the dust is compact (r<19 pc) at 
mid-IR wavelengths. On these two counts at least, it satisfies the requirements of the AGN 
unified model. On the other hand, our modeling shows that the accretion-powered nucleus 
of NGC 1097 is incapable of heating the dust efficiently enough to produce the observed 
mid-IR emission. The torus demanded by the unified model, if it exists at all, should be 
about a factor of two weaker in the mid-IR than the entity that we have detected, and 
have a different SED. The dust that we observe is not heated by the central engine and is 
therefore not necessarily d irectly associated with the AGN, consistent with the prediction of 
Elitzur fc Shlosmanl (120061 ) that tori should disappear at L < 10 42 erg s _1 . 



This does not entirely rule out, however, that NGC 1097 hosts a dusty torus. If a cool, 
IR-faint torus previously existed, Fig. 0] shows that the young nuclear star cluster is capable 
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of heating the dust so that it emits quite strongly in the mid-IR. Alternatively, the starburst 
itself could be the obscurer of the unified model. For this to be plausible, the dust associated 
with the starburst must provide only partial covering and star formation must be common in 
the innermost regions of at least low-luminosity AGN. The geom etry of the dust is n ot well 
constrained by the observations and modeling presented here, but iPrieto et al.l (120051 ) detect 
cold, dusty filaments spiraling down into the central 10 pc of N GC 1097, apparently in a 
thin disk geometry, and a kine matic study by Fathi et al.l ( 20061 ) shows material streaming 
into this central 10 pc region. iFathi et al.l (120061 ) suggest that the infalling material may 
have triggered the nuclear star formation, and disks of fairly recently-form ed stars seem 
quite common in the central ~50 pc of Seyfert galaxies (IDavies et al.ll2006al ). Obscuration 
by dusty starbursts seems at odds with the fact that low-luminosity AGN appear to have 
unobscured nuclei (about as many of the 1 7 LINERs with known nuclear UV point sources 



are type 2s as type Is; iMaoz et al.l 120051 ). but could still apply to some fraction of the 



population. An observational assessment of the frequency of very small-scale star formation 
in these objects would be valuable in this context. 

To the best of our knowledg e, small-scale nuclear dust emission has been searched for in 



only a handful of other LINE Rs (jChary et al.ll2000l : iGrossan et al.ll200ll ; iPerlman et al. 



2001 



Whvsong fc Antonuccill2004l ). of which such emission was detected only in M81. IGrossan et al 
( 120011 ) note that the SED of M81 resembles those of Seyfert galaxies and quasars, suggesting 
that the nuclear 10 /im emission in that object is also from a warm, central dust concentra- 
tion. These results show that at least some LINERs do have dusty nuclei, but it is not yet 
clear whether these galaxies, particularly NGC 1097 with its nuclear starburst and infalling 
material, are unusual cases. High spatial resolution infrared observations of a larger sample 
of objects would help to establish whether or not this is true of AGN-containing LINERs in 
general, elucidate the mechanisms governing the presence or absence of the dust, and clarify 
how these low-luminosity objects may fit into unified schemes of AGN. 



4.2. The lack of strong PAH emission at the nucleus 

As discussed in § 13.21 3.3 /im PAH emission is not detected in the nucleus of NGC 1097, 
although a strong 11.2 /xm band is seen in the larger-aperture IRS spectrum (Figures [5] & 
If)]). Most of the L-band flux in the smaller aperture must arise in material that does not 
produce a 3.3 /im band and whose emission becomes more important in smaller apertures. 
Likely sources of such emission include hot dust around the AGN and nuclear stellar cluster, 
and the central engine itself. The intrinsic AGN/RIAF contribution to the L-band spectrum 
is not taken into account in Fig. HI but it would have to account for most of the total in 



-25 - 



order to "dilute" the PAH emission enough to avoid its detection in the GNIRS spectrum. 
The model fits suggest that this is unlikely (most of the L-band emission is thermal), even 
taking into account the simplicity of the starburst model, but some contribution is certainly 
possible. As for hot dust, that dust must be depleted (relative to the surrounding regions) 
in PAH molecules able to emit at 3.3 /xm. 

Both the AGN and starburst may contribute to suppressing the PAH emission. Calcula- 
tions suggest that PAH mo lecules should be destroyed close to luminous AGN; in particular, 
Siebenmorgen et al.l (12004 ) estimate that the PAH band carriers should be destroyed by soft 
X-ray fluxes > 10 3 erg s _1 cm -2 . Extending the spectral model of iNemmen et al.l (120061 ) and 
correcting for extinction, the intrinsic AGN luminosity from 0.02 to 2 keV is 6.2 x 10 40 erg s^ 1 . 
At a distance of 150 pc from the nucleus — the most distant regions covered by the IRS 
spectrum — the flux is then only about 10 _2 erg s _1 cm -2 , consistent with the appearance of 
the PAH bands in the Spitzer data. Even at only 1 pc, though, the (current) soft X-ray flux 
from the AGN is a factor of a few below the PAH destruction threshold. 

While the starburst does not provide much X-ray emission, it is luminous at UV wave- 
lengths. The 10 6 M(7), 10 6 year-old sta rburst has an intrinsic luminosity above 0.02 keV of 



2.1 x 10 42 erg s 1 (ILeitherer et al.lll999l ). Because of the starburs t's softer radiation field, th e 



average flux must exceed ~ 10 4 erg s 'cm 2 to destroy PAHs (jSiebenmorgen et al.l I2004J ) . 



The starburst would be unlikely to destroy the distant molecules detected in the IRS spec- 
trum, but it could destroy the PAHs located within ~1 pc. 

Observations of other stellar clusters of similar mass and age support the idea that the 
stellar clus ter itself is partly responsible for the lack of PAH emission in its immediate sur- 
roundings. iTacconi-Garman et al.l (120051 ) recently showed that, while the peak 3.3 /zm band 
flux occurs in the general nuclear regions of the starburst galaxies NGC 253 and NGC 1808, 
the feature-to-continuum ratio actually reaches a minimum at the locations of maximum star 
formation activity in those nuclei. This is ascribed to photoionization or photodestruction 
of the PAH molecules in the radiation field of the stellar clusters. The strong PAH em ission 
from the massive, young clusters in galaxies like He 2-10 (lAlonso-Herrero et al.l 120041 ) might 
simply arise in diffuse material outside the immediate vicinity of the star clusters but still 
included in the 1.6" slit used for that wo rk; a similar explan ation has been suggested for 
the PAH bands in the Antennae galaxies (jSnijders et al.l 120061 ) . This would account for the 
strength of the PAH bands in the Spitzer spectrum of NGC 1097 and their weakness at the 
nucleus. 

To summarise, the 3.3 /zm band has b een a useful tracer of star formation on fairly small 
scales ( ~ 60 - few hundred p c) in AGN ( llmanishil |2002| . 120031 ; iRodrfguez-Ardila &: Viegas 



20031 ; llmanishi fc Wadall2004l ). However, if NGC 1097 turns out to be typical, the combined 
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effects of AGN continuum and the radiation field of stellar clusters themselves probably 
mean that it is unlikely that the 3.3 /im feature will be as useful a diagnostic of vigorous 
star formation occurring on the smallest (few parsec) scales. 



5. Conclusions 

We have performed mid-IR imaging and L-band spectroscopy of the central region of 
NGC 1097, at <0.5" angular resolution, to search for thermal dust emission from the active 
nucleus and to test a potential diagnostic of small-scale nuclear star formation. A compact 
nuclear source is detected at both 11.7 and 18.3 /mi, with a luminosity wel l in ex cess of 



that predicted by the accretion disk + RIAF + jet model of iNemmen et al.l (120061 ). This 
mid-IR "excess" is strong evidence for thermal emission from dust in the nucleus of this 
low- luminosity AGN. 



Fitting clumpy obscuration models to this and other 



that the star cluster detected by IStorchi-Bergmann et al.l (120051 ). not the central engine, is 



ligh-r esolution IR data suggests 



the dominant contributor to the heating of the dust in the nucleus. In fact, tori heated solely 
by low-luminosity AGN — if they exist at all — are likely to be relatively faint even in the 
mid-IR close to where their thermal emission peaks. This may mean that the mid-IR emission 
that we detect is not from the torus of AGN unified schemes, in accordance with certain 
predictions that tori should not exist in AGN with luminosities as low as that of NGC 1097. 
Alternatively it is possible, although not required by the data, that dust associated with the 
small-scale nuclear starburst is itself the obscuring material of the unified model. Further 
observations of IR emission and signatures of nuclear star formation in low-luminosity AGN 
would help distinguish between these possibilities. 

Despite the likely energetic importance of the starburst, the 3.3 /im PAH band is not 
detected in the central 19 pc of this galaxy, although strong 11.2 /iin PAH emission is 
observed in the surrounding areas. Destruction and/or ionization of the PAH molecules in 
the radiation field of the stellar cluster itself are likely to account for the non-detection of 
the 3.3 /im band. These factors may limit the use of PAH features as extinction-insensitive 
tracers of very compact nuclear starbursts associated with dusty tori in other low-luminosity 
AGN. 
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